Electronic Raman scattering in quantum dots revisited 
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We present theoretical results concerning inelastic light (Raman) scattering from semiconductor 
quantum dots. The characteristics of each dot state (whether it is a collective or single-particle 
excitation, its multipolarity, and its spin) are determined independently of the Raman spectrum, in 
such a way that common beliefs used for level assignments in experimental spectra can be tested. 
We explore the usefulness of below band gap excitation and an external magnetic field to identify 
charge and spin excited states of a collective or single-particle nature. 
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I. INTRODUCTION 

Raman scattering in semiconductor structures was de- 
vised more than twenty years ago by Burstein et. al. 
as a powerful tool for the identification of electronic 
excitations^. To the best of our knowledge, experiments 
on Raman scattering in quantum dots were performed 
mainly before 1998iSii4 The lack of theoretical calcula- 
tions for the relatively large dots used in the experiments 
(dozens of electrons per dot) made the experimental re- 
sults less conclusive. 

A second handicap of the experiments reported in Refs. 
is related to the fact that they explored the con- 
ceptually difficult resonant regime, where the incident 
photon energy is in resonance with an electronic state in 
the conduction band. In contrast, the existing (qualita- 
tive) theory of Raman scattering is expected to be valid 
only well away from resonance^. It is usually called the 
off resonance approximation (ORA), and is well known 
for missing out the single-particle peaks in the Raman 
spectrum^, which are particularly important in the reso- 
nant regime. In Ref. laser excitation energies 40 meV 
above band gap were used to identify collective states. 
The physics of Raman scattering under these conditions 
is expected to be still more complex because of the sud- 
den increase of level widths with the opening up of the 
channel for spontaneous emission of longitudinal optical 
(LO) phononsi. The positions of collective peaks were 
computed by means of the ORA&, which constitutes a 
nice example of how one can get reasonable results with 
a theory that does not work in this regime. 

In the present paper, we review a set of theoretical 
results on Raman scattering in relatively large quantum 
doti-vliiJ-L!^. which were motivated by the experiments 
described in Refs. fH^. 

We explore the below band gap excitation regime to 
study how the ORA is reached. In addition, it is shown 
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that this regime is ideal for the identification of the collec- 
tive peaks and certain single-particle excitations (SPEs). 
Polarization rules for collective states, following from the 
ORA, are tested. They are unexpectedly shown to work 
also for the SPEs at zero magnetic field, and to break 
down in the presence of a field. Jump rules for the Ra- 
man peak intensities when the laser excitation energy 
approaches the band gap are also demonstrated. 

Raman spectra computed with resonant excitation 
(but with the laser energy below the threshold for the 
creation of LO phonons) show extreme sensitivity to the 
excitation energy, and a reinforcement of single-particle 
peaks. 

We hope that our detailed findings will motivate more 
experimental work on the Raman spectroscopy of quan- 
tum dots, which is regaining interest as a tool for study- 
ing small (with less than 7 electrons) self-assembled 
dotsii 



II. GENERALITIES 
A. The experimentalist point of view 

A typical experimental setup for Raman measurements 
is sketched in Fig. 2] Usually, linear light polarization is 
used, although circular polarization would be very conve- 
nient in certain situations. We take the incident electric 
field to be in the plane of the dot, which is assumed quasi 
twodimensional. 

The backscattering geometry, in which the incident 
and scattered light beams go through the same optical 
fiber (</>i = 4>f), is also very common. 

Two kinds of measurements are usually performed in 
order to separate charge and spin excitations. In the so 
called polarized spectrum, (a), the final and incident elec- 
tric fields are parallel, Ei \\ Ef. Whereas the depolarized 
spectrum is taken with Ei -L Ef, (b). Notice that oblique 
incidence (4>i, (f>f ^ 0) is used to excite multipole final 
states. 




FIG. 1: Geometry of the Raman experiment: (a) polarized, and (b) depolarized geometry. 



B. The theorist point of view 

The theoretical representation of a Raman process 
emerges from the second-order perturbative result for the 
transition amplitude^ 4 , A${\ 

v (f\Hj_ r \int)(int\H^ r \i) 
fl ~ 4f hvi - (E mt - Ek) + iT int ' 1 ' 
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A schematic representation is given in Fig. The 
Raman scattering process can be viewed as the result 
of virtual transitions via all possible intermediate states. 
In the experiments, the incident photon energy, hfi, is 
close to the effective band gap, so that we can restrict 
the sum to intermediate states that contain an additional 
electron- hole pair. 

In the experiments, the temperature, T, is usually be- 
low 2 K. We present calculations at exactly T — K, in 
which only the ground state of the iV-electron quantum 
dot is initially populated. This will be our initial state, 

The intermediate electronic states, \int), as mentioned 
above, are states with an additional electron-hole pair. 
This pair is (virtually) recombined leading to a photon 
with energy hvf. The electronic system ends up in a fi- 
nal state, |/). The energy difference, AEf — Ef — Ei, is 
the Raman shift. Notice the width, Ti nt , entering the ex- 
pression for the transition amplitude in Eq. (JIJ ■ We will 
consider intermediate states in an energy window of 30 
meV above the band gap in order to neglect abrupt vari- 
ations of Ti n t due to LO phonons. Ti nt will be fixed phe- 




FIG. 2: The quantum mechanical transition amplitude and 
its interpretation in terms of virtual transitions. 



nomenologically to a constant value, equal to 0.5 meV, 
which is due mainly to pair recombination. 

From the transition amplitude one computes the cross 
section: 

da~J2\ A fi\ 2 ^( E i + ^i-E f -hu f ). (2) 
/ 

Notice the energy conservation forced by the delta func- 
tion, which implies that the Raman shift is equal to the 
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FIG. 3: The calculation scheme used in the paper. 



photon energy loss. In our calculations, we replace the 
delta function by a Lorentzian: 



where the width of final states will also be assumed con- 
stant, Tf — 0.1 meV. A 30 meV excitation window for 
final states will be used in order to avoid considering 
phonon processes. 



C. Calculation machinery 

The calculation scheme, explicit expressions, and ex- 
planations can be found by the interested reader in pa- 
pers |l0ll5l | . We shall give in this section an overview of 
the conceptual and numerical procedure with the aid of 
Fig. E| 

The dot is modelled by a quantum well with hard 
walls in the z direction (the growth direction of the het- 
erostructure), and a soft harmonic confinement (Hloo = 
12 meV) in the x, y plane. The basis functions used to 
describe one-particle states in the dot are constructed as 
products of harmonic oscillator, infinite well, and spin 
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functions. 

To save time in the calculation of many-electron wave- 
functions, we computed the matrix elements of Coulomb 
interactions, (a, (3\l/r\~/5), where a, /3, 7,<5 are arbitrary 
one-particle states, and stored them in a computer file. 
This calculation takes around 7 days in a personal com- 
puter. The matrix elements are loaded into the com- 
puter at the beginning of a calculation, allowing us to 
solve the nonlinear integro-differential Hartree-Fock (HF) 
equations for 42 electrons in a few minutes, or to compute 
all of the intermediate states entering a Raman process 
(around 10 000) in a few days. The HF equations for 
holes include the electrostatic field created by the back- 
ground electrons in the dot, and take account of valence 
band mixing effects, as described by the Kohn-Luttinger 
Hamiltonian 16 . Up to 6 quantum- well sub-bands are con- 
sidered in the HF equations for holes. 

The final states of the Raman process are excited states 
of the TV-electron system. The intermediate states, on 
the other hand, are states with N + 1 electrons and one 
hole. Both kinds of states are computed by means of 
random-phase-approximation (RPA) like ansatzs for the 
wavefunctionsii, as illustrated in Fig. [3] It was already 
mentioned elsewhere^ that, in our opinion, the main lim- 
itation of using these functions in the present context is 
not related to the well-known lack of correlation effects in 
the RPA, but to the incomplete description of the density 
of energy levels in intermediate and final states. 

Final states are labelled by the quantum numbers AL 
and AS Z , which refer to changes (with respect to the 
ground state) in the total angular momentum and total 
spin projections along the z axis. Borrowing a terminol- 
ogy from Nuclear Physics, we speak about monopole exci- 
tations when AL = 0, dipole excitations when AL = ±1, 
and quadrupole excitations when AL = ±2. States are 
further classified by the degree of collectivity, computed 
with the help of energy-weighted sum rulesit For charge 
monopole excitations, for example, we have: 

T, AE f\ D °M 2 - 2 ^ £<V|A>. (4) 

/ \<fJ,F 

The left hand side of this equation contains only many- 
particle magnitudes, i.e., the excitation energies and the 
matrix elements of the monopole operator between the 
initial and final states, D ^ (its explicit definition can 
be found in Ref. ITfih . The right hand side, however, 
can be evaluated in terms of the occupied HF orbitals, 
A. The magnitude p is the cylindrical coordinate in the 
plane. We will conventionally say that the final state / is 
a collective state if AEf |D"J 2 is greater than 5% of the 
right hand side of Eq. |0J. The same statement applies 
for other multipolc excitations. 

Concerning the spin quantum numbers, the formalism 
leads to ladders of final states characterized by AS Z . In 



order to discriminate between different total spin states 
within a single ladder, we compare the energy differ- 
ences of states in adjacent ladders with the correspond- 
ing Zeeman splitting. We talk about charge excitations 
above the ground state when the total spin does not vary, 
AS* = 0, and spin excitations when AS = ±1. The dis- 
tinction between collective and single-particle states is 
done on the basis of energy-weighted sum rules, which 
can be written for any multipolarity and AS*. Collective 
excited states with AS* = are called charge density ex- 
citations (CDEs), whilst collective states with AS = ±1 
are called spin density excitations (SDEs). The single- 
particle excitations are further distinguished by their 
charge or spin nature. Thus, we will refer to them as 
SPE(C) or SPE(S). 

The intermediate states, on the other hand, are char- 
acterized by the spin S z of the added electron, and the 
total pair angular momentum, T = l e + lh ~ rrij . 

Once the intermediate and final states of a Raman pro- 
cess are computed, the matrix elements of the electron- 
radiation interaction Hamiltonian are evaluated in terms 
of the coefficients entering the ansatz for the many- 
particle wavefunctions and the coefficients of the HF ex- 
pansion. A complete Raman calculation, which requires 
the sum over all the intermediate and final states (previ- 
ously stored into the computer), takes around one day in 
a personal computer. 



D. The OR A 

The ORA is a limit in which the expression for the am- 
plitude of a Raman process is simplified. In this limit, 
only the initial and final state wavefunctions enter the 
expression for A ji . The sum over intermediate states dis- 
appears from it. This means that, within the ORA, one 
can not explain phenomena such as intermediate state 
resonances in Raman processes. 

The explicit derivation of the ORA for Raman scat- 
tering in quantum dots can be found in Ref. 11]. The 
assumptions for its derivation are basically two: (i) The 
laser excitation energy is far away from any intermedi- 
ate state energy, in such a way that we can neglect the 
dependence of the denominator of Eq. on Ei nt , and 
(ii) There is an energy window in the intermediate states 
(from E gap to approximately E gap + 40 meV) where the 
variations of Ti nt can be neglected and the completeness 
relation is practically fulfilled: Ylint Iwt) (int\ ~ 1. Un- 
der these assumptions, one can write: 



A% RA ^ (f\H+_ r H-_ r \i). (5) 

Using the explicit expression for Lf e -r , one arrives at the 
ORA formula: 
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FIG. 4: Calculated Raman spectra in different channels. The incident laser energy is hvi = E gap — 5 meV. 
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A few important conclusions may be derived from Eq. 
©. First, we notice that only collective final states have 
a nonvanishing amplitude in this approximation. The 
SPEs play no role in the Raman process. Indeed, by 



expanding the exponential in Eq. one can obtain 

an alternative expression for A^ t RA in terms of multipolc 
operators^. Only final states having nonzero matrix el- 
ements of multipole operators will contribute to A9^ A . 
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A second important conclusion is related to the spin 
selection rule for Raman scattering. Notice that, in Eq. 
©, multipole operators that do not alter the spin quan- 
tum numbers of the initial state are multiplied by the 
factor Si ■ Ef. This means that peaks corresponding to 
charge operators will appear in the polarized geometry. 
On the other hand, multipole operators that modify the 
spin are multiplied by the factor EiXEj and, consequently, 
Raman peaks corresponding to spin excitations are ex- 
pected to be seen in the depolarized geometry. 



III. RESULTS 

In the next two subsections, we present results for Ra- 
man intensities in quantum dots for laser excitation en- 
ergies below and above band gap. In the former case, 
no experimental measurements have been performed so 
far. We expect weak Raman signals in this regime, but 
there are also many advantages such as, for example, the 
absence of a luminescence background, a smooth depen- 
dence of peak intensities on the excitation energy, etoi£. 
On the other hand, in the above band gap excitation 
regime, our calculations are performed for an excitation 
window ranging from E gap to E gap + 30 meV, where E gap 
is the effective quantum dot band gap. This is sometimes 
called the extreme resonance region. 



A. Raman spectra with below band gap excitation 

The salient features of the Raman spectra in this ex- 
citation regime are summarized below. 
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FIG. 5: \A fi \ 2 for the monopolar CDE and two SPEs(C) 
with excitation energies 17.8 and 18.4 meV, respectively. 



parison. The laser energy is 5 meV below the band gap, 
and the incident (and backscattered) angle is equal to 
20°. Charge monopolar, charge quadrupolar, and spin 
monopolar peaks (with no spin flips) exhibit comparable 
magnitudes. Dipolar (not shown) or spin-flip excitations 
lead to Raman peak intensities one or two orders of mag- 
nitude lower than charge monopolar peaks. 



1. Dominance of monopole peaks 



There are two reasons for the final state monopole ex- 
citations to be the dominant peaks in the Raman spec- 
trum. First, Raman scattering proceeds via the (virtual) 
exchange of two photons. Selection rules dictate that the 
variation of the angular momentum (with respect to the 
initial state) should be preferably AL — 0, ±2, etc. Sec- 
ond, the band-orbital factor in the matrix elements of 
H e _ r (see Fig. [3J provides roughly a factor (qr) 1 every 
time a pair is created or annihilated with a pair angular 
momentum I. As quantum dot dimensions are typically 
r ~ 100 nm^2> and the laser wavelength is A ~ 700 nm, 
pairs with I = dominate the process. 

On the other hand, final state spin excitations in which 
there is a spin flip with respect to the initial state are de- 
pressed, as discussed elsewhereiS. The Raman amplitude 
turns out to be proportional to the minority component 
of the Kohn-Luttinger hole wavefunction. This means 
that monopole spin excitations in which AS Z = should 
be the dominant peaks in the depolarized geometry. 

We show in Fig. 01 some spectra in different angular 
momentum and spin channels for the purpose of com- 



2. Smooth dependence of peak intensities on the excitation 
energy 

The intensities of individual Raman peaks show a 
smooth dependence on hvi when the latter is below the 
band gap. This is simply understood from the Eq. for 
the transition amplitude. From the experimental point of 
view, it is a nice feature. The identification of individual 
peaks could not be easier. 

According to the ORA, as hv{ moves away from the 
band gap, the collective states should start dominat- 
ing the Raman spectrum. In this way, we can identify 
the collective and single-particle excitations (mainly the 
monopolar and quadrupolar ones) by varying hut. We 
show in Fig. 0the amplitude squared, |^4/i| 2 , computed 
in the polarized geometry, for three charge monopolar fi- 
nal states. One of them is the CDE, and the other two 
are SPEs with excitation energies 17.8 and 18.4 meV, 
respectively. The CDE already becomes the dominant 
peak when hvi is around 3 meV below the band gap. 
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FIG. 6: Polarized and depolarized monopolar Raman spectra 
and their comparison with the density of final state SPEs. 
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FIG. 7: Polarized and depolarized monopolar Raman spectra 
at B = and 1 T. 



3. Correlation between the SPE Raman peaks and the 
density of final state energy levels 

It is natural to expect Raman peaks to be located there 
where there is an agglomeration of final states^. In this 
paragraph, we take a step further and compare the po- 
larized Raman spectrum with the density of final state 
SPEs(C), and the depolarized spectrum with the density 
of SPEs(S). This is an attempt to test the spin selection 
rules, derived from the ORA for the collective excitations, 
in the SPEs. 

The comparison is given in Fig. where we show 
results in the monopolar channel. Unexpectedly, the cor- 
relation is high, particularly in the depolarized geometry. 
It is instructive to comment on a commonly held belief 
that the SPEs appear in both polarized and depolarized 
spectra. Of course, it is true. But it would be better 
to say that the SPEs(C) appear mainly in the polarized 
geometry, and the SPEs(S) mainly in the depolarized ge- 
ometry. This statement is valid at zero magnetic field. 



4.. Breakdown of the polarization selection rules in a 
magnetic field 

In a magnetic field, the selection rules deduced from 
the ORA for the collective states are no longer valid at 
excitation energies close to the band gap. The reason 
is the magnetic field dependence of energies and wave- 



functions of intermediate states entering the summation 
of Eq. QJ. The situation is depicted in Fig. \7\ where 
monopolar Raman spectra at B = and 1 T are drawn. 
The excitation energy is hui = E gap — 2.5 meV. 

Let us define the polarization ratio of a sin- 
gle final state, |/), as the ratio of |^4/i| 2 in 
unfavorable and favorable geometries, i.e., r = 
\Afi (unfavorable) | 2 / 1 Afi (favorable) | 2 . By favorable we 
mean the polarized geometry for a charge excitation, and 
the depolarized geometry for a spin excitation. 

The polarization ratios for the collective states at 
B = are 3 X 1CT 4 for the CDE, and 2 x 1(T 7 for the 
SDE. That is, an almost perfect fulfillment of the ORA 
selection rules in spite of the fact that hv{ is only 2.5 meV 
below band gap. At B = 1 T, however, these numbers 
change to 0.03 and 0.17, respectively. 

The fulfillment of the selection rules for the SPEs is not 
so evident in Fig. because these states, with excitation 
energies between 16 and 19 meV, are not as easily dis- 
tinguishable as the collective ones. From the data used 
to draw this figure, we can compute r for each SPE. The 
average polarization ratios at B = arc 3 x 10~ 3 for the 
SPEs(C), and 2 x 10~ 4 for the SPEs(S). At B = 1 T, 
these numbers become 0.26 and 1.1, respectively. That 
is, the intensities in both geometries become comparable. 
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FIG. 8: Polarized and depolarized monopolar Raman spectra 
at B — 4.5 T as hvi approaches the band gap. 



5. Jump rule at the band gap 

Near the band gap, the intensities of Raman peaks fol- 
low an interesting behavior. Let us focus on the collec- 
tive excitations and compare the intensities of the CDE 
and SDE peaks as a function of hvi. The results for the 
monopolar spectra at B — 4.5 T are shown in Fig. |SJ 

We observe a smooth increase of peak intensity in the 
favorable geometry as h,Vi approaches the band gap. In 
the unfavorable geometry, however, the intensity remains 
practically constant up to the moment when hv% reaches 
the band gap, where there is a sudden variation. We 
call this phenomenon the "jump rule" . SPEs also largely 
follow this rule. 



B. Raman spectra under resonant excitation 

We return now to the B = case, and increase the 
laser energy to values above the band gap. As mentioned 
above, we restrict hvi to the interval (E gapi E gap + 30 
meV) in order to avoid considering phonon effects. 

Many of the properties discussed in the preceding sec- 
tion are still valid in the present context. For example, 
the monopolar and quadrupolar peaks are the most im- 
portant peaks in the Raman spectrum, and the density 
of final state energy levels is correlated with the positions 
of the principal Raman lines. 
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FIG. 9: Monopolar Raman spectra above band gap in the 
polarized geometry. 



In addition, there are certain new features, which char- 
acterize the resonant excitation regime. 



1. Dominance of SPEs 

Raman peaks related to SPEs experience a noticeable 
increase of intensity under resonant excitation. Qualita- 
tively speaking, one can say that Raman scattering pro- 
ceeds through a single intermediate state (the one exactly 
at resonance), which (virtually) decays indiscriminately 
to the collective final state or to the SPEs. As there is 
a large number of SPEs, and they are packed in groups, 
the intensities of the corresponding peaks can be high. 

This statement is illustrated in Fig. [5J where three 
spectra corresponding to hvi = E gap + 5 meV, E gap + 10 
meV, and E gap + 15 meV, respectively, are shown. In 
addition, we observe a non-monotonous dependence of 
peak intensities on hvi. This is a consequence of the fact 
that, when hv% is varied, a different intermediate state 
comes into resonance. 



2. Correlation between the Raman intensities of individual 
states and the density of intermediate state energy levels 

Here, we show how by monitoring the intensities of in- 
dividual peaks as a function of hvi one can obtain infor- 
mation about the density of energy levels in intermediate 
states, at least at low excitation energies above the band 
gap. 
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FIG. 10: Upper panel: \A fi \ 2 for the same three charge 
monopolar states used in Fig. and hvi above band gap. 
Lower panel: Histogram of intermediate state energy levels. 



In Fig. w e follow the same final states used in Fig. 
El and compare the corresponding |A/i| 2 for hvi in the 
interval (E gap , E gap + 30 meV) with the density of energy 
levels. 

The first peak of |^4/i| 2 near 5 meV above the band 
gap signals the beginning of a group of energy levels. A 
second structure is seen near 10 meV above the band 
gap, where there is also a threshold. At higher excita- 
tion energies there is still some correlation between peaks 
in |A/i| 2 and peaks in the density of intermediate state 
energy levels. 
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FIG. 11: |^4/i| 2 for the spin monopolar state with excitation 
energy AEf — 18 meV and the contribution of each individual 
intermediate state to the sum. 



sum. No strong cancellation effects nor strong coopera- 
tion effects are apparent. That is, interference effects in 
Afi are weak under resonant excitation. 



4- Selection rules in a magnetic field 

The polarization selection rules that follow from the 
ORA may be tested for excitation energies above the 
band gap. Surprisingly, they shown to be very well 
obeyed at zero magnetic field, and break down in a mag- 
netic field. The polarization ratio exhibits a strong de- 
pendence on hut and on the applied field. In quality of 
example, let us consider the monopolc SDE. At zero field, 
r is practically zero (1CP 7 when hvi — E gap + 5 meV). 
But at B = 1 T, r varies from 1.1 at the band gap to 0.1 



for hv; = E„ 



2.5 meV. 



We now evaluate the contribution of the different in- 
termediate states entering the summation of Eq. Q in 
resonant Raman scattering. The postulate is that the 
qualitative picture sketched in Sec. IIII B II is correct: 
only those intermediate states whose energies are very 
close to hvi contribute to the Raman amplitude. To ver- 
ify this, we performed different calculations of Afi, re- 
stricting the sum over intermediate states to the window 
{hvi — SE,hi>i + SE). There is qualitatively no change 
in the spectrum when SE is reduced from 10 down to 2 
meV. 

We compare in Fig. ^2 the magnitude |^4/i| 2 for the 
spin monopolar peak with AEf = 18 meV with the in- 
dividual contributions of each intermediate state to the 



IV. CONCLUSIONS 

These theoretical calculations of the Raman spectrum 
of a many-electron quantum dot in zero and non zero ap- 
plied magnetic field have revealed intriguing new features 
and a great sensitivity to the excitation energy. 

For excitation energies below the band gap, we found 
Raman spectra dominated by the collective excitations 
already when hvi = E gap — 5 meV. The peak intensities 
depend smoothly on the excitation energy. In zero field, 
the simple-minded polarization selection rules - CDE po- 
larized/SDE depolarized - are very well obeyed even by 
the SPEs. This means that one can obtain information 
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about the density of final state SPEs from the Raman 
spectra. In an applied magnetic field, these basic selec- 
tion rules break down, particularly when hvi is increased 
up to the band gap, where the Raman intensities of the 
collective excitations in their "forbidden" geometries sud- 
denly increase to appreciable values. This we termed the 
Raman intensity "jump rule" , and may help to determine 
the charge or spin character of a single excitation. 

A quite different situation is encountered for excitation 
at energies above the band gap. The Raman peak in- 
tensities fluctuate considerably as the excitation energy 
sweeps through resonance with intermediate states. In 
this regime, the SPEs dominate the Raman spectrum, 
the Raman intensities of individual peaks correlate well 



with the density of intermediate state energy levels, and 
there is an absence of interference effects. 

It would be interesting now to explore experimentally 
these new features of the electronic Raman spectrum of 
many-electron quantum dots. 
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